Universal Non-Gaussian Velocity Distribution in Violent Gravitational Processes 
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We study the velocity distribution in spherical coUapses and cluster-pair collisions by use of N- 
body simulations. Reflecting the violent gravitational processes, the velocity distribution of the 
resultant quasi-stationary state generally becomes non-Gaussian. Through the strong mixing of the 
violent process, there appears a universal non-Gaussian velocity distribution, which is a democratic 
(equal- weighted) superposition of many Gaussian distributions (DT distribution). This is deeply 
related with the local virial equilibrium and the linear mass-temperature relation which characterize 
the system. We show the robustness of this distribution function against various initial conditions 
which leads to the violent gravitational process. The DT distribution has a positive correlation 
with the energy fluctuation of the system. On the other hand, the coherent motion such as the 
radial motion in the spherical collapse and the rotation with the angular momentum suppress the 
appearance of the DT distribution. 
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I. INTRODUCTION 

Galaxies and clusters of galaxies are typical structures 
formed through the gravity of their own. We would like 
to understand the history and universal characterization 
of these self-gravitating structures. Especially, we are in- 
terested in how extensively the formation process of them 
are involved in the resultant universal structure through 
their gravitational interaction. In these structure for- 
mation of self-gravitating systems (SGS), a cold collapse 
and a cluster-pair collision would be the fundamental pro- 
cesses. Therefore in this paper, we would like to focus on 
such fundamental dynamics disregarding the other non- 
gravity factors. 

The cold collapsing process has been extensively stud- 
ied as a crucial relaxation process of the coUisionless sys- 
tems such as elliptical galaxies where the stellar encoun- 
ters are unimportant. After the violent cold collapse, 
a steady state is generally formed, whose density pro- 
file is well described by the de Vaucouleurs's r^/^ law 
[H Qi IM IJ; 19 ' The spherical averaged density profile is 
found to be p ex r^"' in spherical cold collapses ^.^ij and is 
found to be p ex r^^ — r^^ in the cluster-pair collision [3|. 
The energy distribution is further discussed in t8| . On the 
other hand, the velocity distribution has not yet been ex- 
tensively studied except for the anisotropy of the velocity 
dispersion I1C|. Recently, Merrall et al. have numerically 
showed that the radial velocity distribution in the cen- 
tral region of the bound particles becomes Gaussian^. 
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After the spherical collapse and the cluster-pair collision, 
Kanaeda and Morikawa showed that the velocity distri- 
bution of all bound particles becomes non-Gaussianj^ 
which is described by the superposed-Gaussian. 

From a viewpoint of statistical mechanics, we can- 
not naively expect the ordinary Gaussian distributions 
for SGS, because the the long-range interaction of 
gravity apparently violates additivity, which is the ba- 
sic standpoint of the ordinary Boltzmann statistical 
mechanics 11]. Then the question is whether we can ex- 
pect any universal distribution function for SGS instead 
of the Gaussian distribution. 

As one of the possible explanations for these non- 
Gaussian distributions in the stationary state with large 
fluctuations of intensive quantities such as temperature. 
Beck and CohenT2| proposed the superstatistics. Ac- 
cording to this proposal, statistical properties of the tem- 
perature fluctuations determine overall non-Gaussian dis- 
tributions. A special choice of the fluctuation leads to the 
Tsallis statistics 13J. 

However, it is clear that the non-Gaussian properties 
of SGS are not always observed everywhere in the Uni- 
verse. Moreover non-Gaussian properties of SGS are 
quite diverse in general and we cannot expect the com- 
pletely universal properties of SGS. For example, we 
studied non-Gaussian properties in the self-gravitating 
ring model|l4|. where many self-gravitating particles are 
constrained to move on a circular ring, which is fixed in 
the three-dimensional space. Only at the intermediate 
energy scale where the specific heat becomes negative, 
only the halo particles which belong to the intermedi- 
ate energy scale, have shown non-Gaussian and power 
law velocity distribution: f{v) ex v~^. In this model, 
the existence of the halo particles plays an essential role 
in the appearance of non-Gaussian distributions. In our 
present paper therefore, we would like to clarify in which 



aspect of SGS and under which conditions the universal 
non-Gaussian properties are observed. 

Non-Gaussian properties would be more naturally ex- 
pected for the coUisionless stage of the evolution (t < trei ) 
than the later coUisional stage of SGS (t ^ tr-ez), where 
trei is a local two-body relaxation timejl^ (In this paper, 
trei is defined by Eq.(|2Jl). This is because in the former 
stage, the system is not thermally relaxed and the local 
equilibrium is not yet established. A long range dynami- 
cal correlation among the whole system develops though 
the long range interaction. This property would yield 
non-additivity of the system and manifest deviation from 
the ordinary Gaussian distributions. On the other hand, 
in the later stage, the local equilibrium is established 
through two-body encounters and the situation is similar 
to the ordinary statistical mechanics which admits Gaus- 
sian distributions. Thus we would like to concentrate on 
the coUisionless stage of SGS in this paper. 

In section ^ we begin our study with typical simu- 
lations for spherical cold collapses and cluster-pair colli- 
sions. We find the same form of non-Gaussian velocity 
distribution in both cases, and then we explore four dif- 
ferent models of the superposed-Gaussian distributions 
to describe this velocity distribution. In section IIIII 
by analyzing the numerical data we show that the non- 
Gaussian velocity distribution observed in our simulation 
are well described by the "Democratic Temperature(DT) 
distribution" . This DT distribution is consistent with the 
fact that we observe the linear relation between the tem- 
perature and the inner mass. In section llVl we study the 
universality of this DT distribution and show that the 
mixing property under the violent gravitational process 
is the essence for the appearance of the DT distribution. 
The effect of coherent motion in the velocity distribution 
is discussed in section[V] The last section lVIl is devoted to 
the discussions and further developments of the present 
work. 



II. THE VELOCITY DISTRIBUTIONS IN 
N-BODY SIMULATION 



A self-gravitating N-body system is described by the 
following Hamiltonian: 



can be written in the following form: 
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where r^ and pi are respectively the position and the 
momentum of the ith particle, m is the mass of the each 
particle, and G is the gravitational constant. For numer- 
ical simulations, we have to introduce a cutoff parameter 
e. For the units of length, mass, and time, we use the 
initial system size R, the total mass M :— Nm, and 
the initial free fall time i// := y/R^/{GM), respectively. 
The local two-body relaxation time trd in this simulation 






0.065cr3(r) 
G2p(r)m ln(l/e)' 
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where a(r) is a velocity dispersion and p(r) is a mass 
density. We use a leap-frog symplectic integrator on 
GRAPE-5, a special-purpose computer designed to accel- 
erate N-body simulations [16J. In all runs, the numerical 
errors in the total energy |Aii^/£'o| have been controlled 
to be less than 10~^. 

In this section, we focus on the velocity distribution in 
the violent gravitational process by N-body simulation. 
Especially we consider two fundamental processes of vio- 
lent gravitational dynamics; a spherical cold collapse and 
a cluster-pair collision. 



A. Spherical Cold collapse process 

We first glance at a typical example of a spherical cold 
collapse process (run SC in Tabled). All particles are 
homogeneously distributed within a sphere of radius R. 
The system is composed of iV = 5000 particles and we 
set the vanishing virial ratio \2K/W\ — initially, where 
K and W are the kinetic energy and the potential energy 
of the whole system, respectively. 

Figd shows snapshots of particle distributions at dif- 
ferent times. Most particles rapidly collapse into the cen- 
ter within the free-fall time t ~ i//. After this collapse, 
some particles obtain positive energy and escape from the 
system, and the rest particles remain bounded and grad- 
ually expand leaving a tight core at the center. This is a 
typical process of the formation of the core- halo structure 
for SGS. 

We now focus on the velocity distributions of the parti- 
cles. Since we are interested in the global property, we ex- 
tract the one-dimensional velocity distribution function 
combining all-directional components of velocity distri- 
butions. We will examine the anisotropy in the velocity 
distribution in Sec^l The velocity distribution thus ob- 
tained is shown for bound particles, which have negative 
energy, in Fig|21 This velocity distribution is apparently 
different from the ordinary Gaussian distribution, espe- 
cially in the small velocity region. There is an apparent 
cusp at the center. This velocity distribution has been 
stably observed just after the collapse during the whole 
simulation time (t ~ tff -^ lOOOt//). We emphasize that 
this non-Gaussian velocity distribution is quite universal 
and robust in various cold collapsing process, as we will 
demonstrate briefly. 



B. Cluster- pair collision process 

We turn our attention to another case of the violent 
gravitational process: The cluster-pair collision. Because 
we would like to clarify the basic process, we choose the 
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FIG. 1: A snapshot of a spherical cold collapse case (run SC 
in Table EJ: top left (t = 0), top right (i = li//), bottom left 
{t — 5tff), and bottom right (i = lOOi//). 
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FIG. 2: A velocity distribution of a spherical cold collapse 
(run SC in Table HJ at t = lOi//. Three best-fitted Gaussian 
distributions are superposed. Gauss (all) is fitted by use of 
the all data. Gauss (small) and Gauss (large) are fitted by 
use of data included for only small velocity and only large 
velocity, respectively. 



simplest head-on collision (run CC in Table llT|) . Each 
cluster has an equal number of particles and all particles 
are homogeneously distributed in each sphere of radius R. 
The initial velocity distribution is set to be Gaussian and 
the initial virial ratio is 1. We set the initial separation 
of the pair to 6R along the x axis. 
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FIG. 3: A snapshot of a cluster collision model (run CC in 
Table |n|l: top left (f = 0), top right {t = 5t//), bottom left 
{t = 20t//), and bottom right {t = 500t//). 



FigO shows the snapshots of particle distributions for 
the run CC at different times. The cluster-pair collides 
at time t ~ 20i// and then gradually merges into a single 
cluster. The velocity distribution profile after this merg- 
ing is shown in Fig0] The velocity distribution profile 
is very similar to the previous spherical collapse case(run 
SC) . The excess of velocity distribution at the small value 
is prominent as before. 

In both cases, we notice that the temperature or the 
velocity dispersion is different from place to place; the 
center of the core is much hotter than the outskirts of the 
system. Therefore it is clear from the beginning that the 
obtained velocity distribution cannot be fitted by a sin- 
gle Gaussian distribution. Thus we propose to describe 
these systems by the superposition of Gaussian distribu- 
tions with various temperatures. This consideration is 
very natural because the coUisionless SGS allows the co- 
existence of many temperatures. This is because (a) the 
coUisionless SGS is not yet thermally relaxed, and (b) 
the virial relation implies that the system shows negative 
specific heat, which spontaneously yields inhomogeneous 
temperature. 

In general, superposition of Gaussian velocity distribu- 
tion fsci'^i) has the form: 
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where A is a normalization constant. The parameter 
function M can be any function of temperature. How- 
ever now in the present spherical symmetric case, we can 
choose it to be the inner mass: Mr := /q p (r') Airr'^dr' . 
This choice is not at all trivial and will be demonstrated 
briefly. Any special choice of the weight function for 
the superposition yields in general arbitrary distribution 
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FIG. 4: Same as Fig|21 but for a cluster-pair collision model 
(run CC in Table IHll at t = 500t//. 



functions. In order to avoid such arbitrariness, we will 
choose the simplest weight functions. Here we consider 
four natural models of the superposed-Gaussian distribu- 
tion. 

• Model 1 (GT) Gaussian-weighted superposition of 
Gaussian distributions with various temperatures. 

The temperature distribution is the Gaussian with 
the dispersion t, and GT distribution fGT{vi,T) 
takes the following form: 
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where Erf(a:) is the error function: 

Erf(2:) := —^ \ dte'^^ 
Vt^ Jo 

• Model 3 (Gcr ) Gaussian-weighted superposition of 
Gaussian distributions with various standard devi- 
ations a := Vt. 

The distribution of the standard deviation is the 
Gaussian with the dispersion 7^, and Gcr distribu- 
tion fGij{viil) takes the following form: 
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where r(a;) is a Gamma function and 
pFg(ai, . . . ,ap; f3i, . . . , (3q; z) is the generalized 
Hypergeometric function. 

• Model 2 (DT) Equal-weighted superposition of 
Gaussian distributions with various temperatures. 

This case corresponds to dT / dM — const, in terms 
of the parameter function M, and we call the ve- 
locity distribution in this case as "democratic tem- 
perature distribution" /^(^^(wi): 
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where Kq{z) is the modified Bessel function. 

• Model 4 (Do- ) Equal-weighted superposition of 
Gaussian distributions with various standard de- 
viation a. 

Since da/dM — const., Da distribution foaivi) 
has the following form: 



o Jo V27rcr' 

^ n0,vf/{2a^)) 

2V2^a ' ^ ' 

where T(z,p) is the incomplete Gamma function. 

In the next section, we examine the above four models 
in relation with our non-Gaussian velocity distributions. 



III. DEMOCRATIC TEMPERATURE (DT) 
DISTRIBUTION 



As shown in the preceding section, the velocity distri- 
bution functions for both processes, spherical cold col- 
lapses and cluster-pair collisions, take the same non- 
Gaussian form. Using the most natural four models of 
superposed-Gaussian distributions, we try to fit our ve- 
locity distribution function and to extract any universal 
character of SGS. 

The fitting results are shown in FigISland|ni Using the 
X^ analysis, the best fit model turns out to be the DT 
distribution in both cases. 
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FIG. 5: A linear-log plot for velocity distribution for the case 
of a spherical cold collapse (run SC) at i = lOi//. The 
best fit Gaussian and the four superposed-Gaussian model 
(Eqs. (131,0,®, and Q) with a best fit parameter are super- 
posed. The x^ of each model is 0.00095(Gauss), 0.00037(DT), 
0.00048(GT), 0.0016(Da), and 0.0025(Ga). The best fit 
model is a DT distribution. 
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FIG. 6: Same as Fig|3 but for a cluster-pair collision 
model (run CC) at f = 500i//. The x of each model 
is 0.0030(Gauss), 0.0021(DT), 0.0030(GT), O.OlO(Dcr), and 
0.015(Gcr). The best fit model is a DT distribution. 
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FIG. 7: A linear-log plot for a velocity distribution of each 
shell for the case of a spherical cold collapse (run SC). We 
divided the whole bound particles into 10 shells including an 
equal mass and calculated the velocity distribution by use of 
the data from t = lOt// and t — 20i//. 
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FIG. 8: Same as Fig[7| but for a cluster-pair collision (run 
CC). We calculated the velocity distribution by use of the 
data from i = 500i// and t — 510t//. 



In order to investigate the structure of velocity space 
for bound particles in detail, we divide the whole par- 
ticles into several shells with equal number of parti- 
cles. We introduce the inner mass coordinate Mr :— 
47r/p dr'r''^p{r') and consider the averaged quantities 
within each shell as local variables. 

The velocity distribution of each shell shows almost 



Gaussian (Figd and |H1) . These figures show that each 
shell has a different temperature and suggest that the 
previous non-Gaussian distribution can be described by 
the superposition of Gaussians with various tempera- 
tures. The temperature is the highest at the most inner 
shell and monotonically decreases toward outside shells. 
The local velocity dispersion < u^ > is plotted against 
the inner mass Mr in Fig|^ and ^| In both cases, it 
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FIG. 9: A mass dependence of velocity dispersion (< 
v'^{Mr) >) for the case of a spherical cold collapse (run SC). 
The mass is normalized by the mass of whole bound particles 
Mb- The velocity dispersion is time averaged and the error 
bar is calculated from the time fluctuation. 
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FIG. 10: Same as FigEl but for a cluster-pair collision (run 
CC). 



is remarkable that the velocity dispersion decreases lin- 
early in the inner mass Mr. We emphasize here that this 
linear relation is quite robust and universal. Actually 
this relation is observed in our two processes with wide 
range class of initial conditions, which will be further dis- 
cussed in the subsequent sections. If we carefully look at 
FigEl the temperature gradient in the inner region grad- 
ually reduces in time, while it remains constant in the 
outer region. Since the mass density at the inner region 
increases in time and eventually the effect of two-body 
relaxation turns on, such a tilt of the temperature gradi- 
ent is thought to be caused by this coUisional effect. A 
careful check of our numerical simulation supports this 
interpretation. Therefore the linearity of the temperature 
gradient is thought to be the property of the coUisionless 
SGS. Further collisional effect would eventually yield the 
uniform temperature distributions. 

The almost Gaussian velocity distribution in each 



shell guarantees the form of Eq.Q and the linearity 
of the temperature-mass relation leads to the relation 
\dM/dT\ — const, in Eq.©, which indicates the DT dis- 
tribution for the velocity distributions. This is perfectly 
consistent with the result that the DT distribution is the 
best fit to the velocity distribution among several models 
which we examined. 

Besides this linearity in the diagram of the local veloc- 
ity dispersion < v^ > and the inner mass M^, we observe 
another prominent fact that the virial relation between 
the potential energy and kinetic energy holds locally at 
each spatial region. Using the local kinetic energy Kr 
and the local potential energy Wr inside the radius r, we 
can locally define the virial ratio \2Kr/Wr\. The time 
evolution of this local virial ratio is depicted in Figllll 
and ^] The value of the local virial ratio stays almost 
unity within the error ±10% everywhere anytime. This 
suggests that not only the whole system but also each 
shell is locally virialized. This fact may be deeply re- 
lated with the robustness of the DT distribution, and 
will be separately discussed in [l^ . 




FIG. 11: The local virial relation for run SC. The ratio 
\2Kr/Wr\ is plotted as a function of Af,.- The virial ratios 
of each shell at different times (i = 5i//, lOi//, 50i//, and 
lOOtff) are superposed. 

In order to examine the stability of DT distribution, 
we calculate the time evolution of the x^ value for DT 
{Xdt) ^^'^ Gaussian distributions (xGauss)' which is de- 
picted in Figllll and 1141 Immediately after the collapse 
or the merging, the velocity distribution always becomes 
the DT distribution during the whole period of our sim- 
ulations. This robustness makes DT distribution one of 
the most characteristic properties of coUisionless SGS. 



IV. UNIVERSALITY OF DT DISTRIBUTION - 
DEGREE OF VIOLENCE 

As we have seen in the preceding section, DT distri- 
bution, supported from the linear < v"^ > - M^ relation, 
is one of the most relevant characteristics of SGS such 
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FIG. 12; Same as Fie llll but for a cluster-pair collision (run 
CC). The virial ratios of each shell at different times (i = 
lOOi//, 200i//, 500t//, and lOOOt//) are superposed. 
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FIG. 13; Time evolution of Xdt of DT distribution and 
"Xxiauss of Gaussian for the case of a spherical cold collapse 



(run SC). 



as the spherical cold collapses and cluster-pair collisions. 
We would like to further clarify the condition for the ap- 
pearance of the DT distribution function. From the be- 
ginning of our numerical experiments, we have implicitly 
chosen the violent processes such as cold collapse and the 
quiet cluster-pair collision. We had in mind the expecta- 
tion that the strong mixing property associated with the 
violent gravitational process would yield the most promi- 
nent characteristics for the collisionless SGS. Actually it 
is well known that a violent mixing becomes an important 
factor to realize the quasi-equilibrium virialized statejlTJ. 

In this section, we will quantitatively demonstrate this 
expectation. Especially we would like to explore the cor- 
relation between the DT distribution and the degree of 
violence of the processes. 

We control the degree of violence of the process by the 
choice of initial conditions in our numerical experiments. 
The initial conditions of all runs are listed in Tabled for 
spherical collapses and Table ITU for cluster-pair collisions. 
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FIG. 14; Same as Fig ll8l but for a cluster-pair collision (run 
CC). 



TABLE I; Initial condition for spherical collapse case. A^ is a 
number of particles and \2K/W\ is an initial virial ratio and 
p oc r~" is an initial density profile and e is a cutoff parameter 
in Eq.lQ. 
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"The kinetic energy is contributed by only the rigid rotation 
around the z axis. 



As the indicator for the degree of violence of the pro- 
cess, we use the total energy fluctuation of particles; more 
violent the process the more rapid the energy mixing of 
each particles. More precisely, the energy fluctuation of 
the ith particle is defined as 



ef(t)-e,(i) 



(8) 



where e^ {t) is the energy of ith particle at time t and * is 
the time averaged value of * during the whole simulation 
period. Then we define the total energy fiuctuation of 



TABLE II: Initial condition for cluster collision case. A'' is the 
total number of particles and the number of particles in each 
cluster is N/2. \2K/W\ is an initial virial ratio of each cluster 
where the kinetic energy K is contributed by only the random 
motion. Krot is a kinetic energy by only the rotation with an 



orbital angular moment around the z axis ( 


Lz) initially. 


run N 


\2K/W\ 


\2Krot/W\{L,) 


CC 5000 
CCLl 5000 
CCL2 5000 


1 
1 
1 



0.1 
0.2 



the system as the sum of (Jeyi for all particles as 
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%/Zb2 :^^CT^., 



(9) 



j=i 



In order to check the efficiency of VAE^ as the indi- 
cator of the degree of violence, we calculate the correla- 
tions with other physical quantities which are naturally 
expected to be connected with the strength of the mixing. 
When violent mixing in the phase space occurs, some of 
the particles obtain enough energy to escape from the 
system through the energy exchange by the potential os- 
cillation. Then the rest of the particles becomes more 
tightly bounded due to the extraction of energy by the 
escaping particles. Therefore both (a) the system size 
and (b) the total energy of the bound particles, Eb{< 0), 
will become smaller after the collapse, provided the mix- 
ing is strong and effective. As the quantitative indicator 
for (a), i.e., the size of the system after the collapse, 
we introduce the half-mass radius divided by it's initial 
value: Rhm- Further, we need to take ensemble average 

< * > of all those quantities for the regular continuous 
indicator for the process. 

As in shown in Fig^l < vAE^ > is apparently cor- 
related with both < Rhm > and < Ei, >. This fact 
supports that the total energy fluctuation < VaW > is 
actually an effective indicator of the degree of strength 
of the mixing. 

Keeping in mind that < y/AE^ > is a good indicator 
of the violence of the process, we now examine the corre- 
lation between the degree of violence in the initial stage 
and the occurrence of the DT distribution for bound par- 
ticles in the quasi-equilibrium stage. We calculated x^ 
values for the fitting of velocity distributions with DT 
distribution and Gaussian distribution for all runs. The 
result is shown in Fig^J where the ratio of Xdt '^^ ^^ 
distribution to Xcauss of Gaussian distribution is plotted 
against < vAEP > for all initial conditions. The value is 
grea ter th an one for the simulations with small value of 

< VAE^ >, which means that the velocity distribution 
is better fitted by Gaussian rather than DT dist ributi on. 
On the other hand, for the larger values of < vAEP >, 
the data is better fitted by DT distribution rather than 
Gaussian. Thus DT velocity distribution is always favor- 
able in the simulation with violent processes. 
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FIG. 15: The half-mass radius normalized by it's initial value 
at t = 0, Rhm and the absolute value of t he tot al energy of 
bound particles, \Ei,\ are plotted against VAEP for all cases 
in Table 13 and Table [H] We take the ensemble average for 
each data with several times. It ranges from t = Wtff to 
t = loot// with Wtff intervals in the spherical case, and 
from t = lOOtff and t = lOQOtff with lOOt// intervals in 
cluster collisions. The error bar is calculated from the ensem- 
ble average. 
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FIG. 16: XDr/xhauss for all cases in TableEland TablelHl x^ 
is a time averaged from t = Wtff and t — IQOi// at Wtff 
intervals for the spherical collapse case, from t = lOOi// and 
t = WOOtff at WOtff intervals for the cluster collision case. 
The error bar is calculated from the time fluctuation. 



This point is further clarified in Fig II 71 where we plot- 
ted the value of Xdt against < vAE^ >. It is appar- 
ent that the Xdt nionotonically decreases for increasing 
vA£^. These results manifestly support our expecta- 
tion that the DT velocity distribution is associated with 
the strength of the violent mixing. 

From these results, we can claim that the DT velocity 
distribution function is universal for various initial con- 
ditions provided that the violent gravitational process is 
included. Moreover, the DT distribution becomes more 
prominent for stronger violent processes. 
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FIG. 17: Same as FigCni but for Xdt- 



V. UNIVERSALITY OF DT DISTRIBUTION - 
DEGREE OF COHERENCE 

We have so far studied the universahty of DT distri- 
bution and its correlation with the degree of violence of 
the dynamical process. Here in this section, we further 
discuss another aspect of universality of DT distribution 
function. 

In the above argument, we have neglected possible 
anisotropy in the velocity space; we have extracted one 
dimensional velocity data set by simply combining all the 
velocity components. However actually in the processes 
of spherical collapse and cluster-pair collision, the veloc- 
ity space is anisotropic in general. Further the anisotropy 
in the rotating collapse (runs SCSl-3 in Tabled is appar- 
ent. 

For example, we first consider the spherical cold col- 
lapse (run SC). FiglTHl shows the measure of validity of 
DT distribution, i.e., the time evolution of Xdt divided 
by xhauss- The radial velocity distribution is better fit- 
ted by Gaussian for a while until t ^ 50t// and after that 
time DT distribution becomes better. On the other hand, 
the tangential velocity distribution is always well fitted 
by DT distribution during all our simulation period. 

This behavior will be related with the coherent mo- 
tion of the particles. More precisely, in the early stage 
of the evolution, there exists a coherent radial motion 
such as a collapse and a subsequent bounce, which even- 
tually decays. This coherent motion does not affect the 
tangential velocity distribution. Moreover this coherent 
motion may inherit some amounts of randomness in the 
initial distribution of particles. Therefore in the early 
stage, only the radial velocity component is dominated 
by the coherent motion with this initial spatially ran- 
dom distribution, and eventually recovers the intrinsic 
DT distribution when the coherent motion decays. 

The above consideration is also applicable for the 
anisotropy in the rotating collapse (runs SCSl-3 in Ta- 
bleHJ. Fig 1 191 shows the measure of validity of DT distri- 
bution, i.e., the time evolution of Xdt divided by xiauss 



in this case. Only the i^x-component shows DT distri- 
bution all the time, and the Uxjt'y-components, the di- 
rection of coherent rotation, never show DT distribution 
during our simulation time. 

Therefore the condition for the appearance of the uni- 
versal DT distribution would be, besides the violence 
of the system, the decay of the coherent motion, which 
keeps the information of initial conditions. 
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FIG. 18: Time evolution of Xdt of DT distribution divided 
by the Xcauss fo^' the case of a spherical cold collapse (run 
SC). 
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FIG. 19: Same as Fie llSI but for the case of a spherical col- 
lapse with the angular momentum L^ (run SCS2). 



VI. CONCLUSIONS AND DISCUSSIONS 

We have studied the velocity distribution function of 
self-gravitating systems (SGS) which experienced vio- 
lent gravitational processes by use of N-body simulation 
method. Especially we have chosen two fundamental pro- 
cesses which mimic the galaxy formation dynamics; the 
spherical cold collapse and the cluster-pair collision. 

In both cases, after the collapses or the collisions, all 
the bound particles form a single stationary state, which 
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are characterized by the local virial relation and the lin- 
earity in the mass-temperature relation. In this station- 
ary state, the velocity distribution is well described by the 
democratic (—equally weighted) superposition of Gaus- 
sian distributions of various temperatures (DT distribu- 
tion). 

This DT velocity distribution is robust against various 
initial conditions such as the change of the particle num- 
ber, the virial ratio, and the density profile. Moreover, 
using the half-mass radius and the root mean square of 
energy variation as a measure of the degree of violent 
mixing, we found a firm positive correlation between the 
DT distribution and the degree of violent mixing. This 
correlation suggests that the DT distribution originates 
from a violent mixing in SGS. 

Carefully examining the anisotropy in the velocity dis- 
tributions, we have also found another condition for the 
appearance of DT distributions. // any systematic coher- 
ent motion exists then it could inherit special initial con- 
dition and could suppress any intrinsic properties such as 
DT distribution for SGS. This coherent motion actually 
exists in the radial velocity distribution in the spherical 
collapse processes, and in the rotational velocity distri- 
bution in the processes with non-vanishing angular mo- 
mentum. 

As a conclusion, we postulate that both the local virial 
relation and DT velocity distribution, associated with the 
linear temperature-mass relation, are universal properties 
of SGS which undergo violent gravitational mixing. 

The origin of the above universal properties should be 
further investigated. One possible candidate would be a 
steady heat flow from the center toward outskirts of the 
spherical system. In our calculation, the two-body relax- 
ation process is initiated at the very center of the core 
region, within our numerical calculation period. Release 
of the gravitational energy at the center would yield a 
steady heat flow toward the outskirts of the system. This 
flow may guarantee the two properties we found in this 
paper and their duration. We hope we can report the 
detail very soon. 

In this paper, we paid attention only to the velocity 



distribution function. However, we need the information 
of matter configuration like density profile to understand 
the full characteristic of gravitationally bound systems 
in quasi-cquilibrium state. We will show that the combi- 
nation of the local virial condition and the linear mass- 
temperature relation leads to a universal density profile 
for gravitationally bound systems in a separate report 
jlSJ I . This information should give a hint for the origin of 
the universality profile of dark matter or elliptical galax- 
ies like de Vaucouleurs law [l|. 

In our previous workp^. we have studied the self- 
gravitating ring model and have obtained a non-Gaussian 
velocity distribution. This velocity distribution shows a 
power law and is apparently different from the results 
obtained in this paper. We believe that this discrepancy 
originates from the difference of the boundary conditions. 
In the self-gravitating ring model, the configuration space 
is compact. On the other hand, in the present simulation, 
the configuration space is open. We would like to exten- 
sively consider the effect of boundary conditions soon. 

In this work, we restricted our models only to the 
spherical collapses and the cluster-pair collisions. How- 
ever actual galaxies may undergo much more diverse evo- 
lution processes including chemical evolution, the envi- 
ronment effects, and a possible special nature of the cold 
dark matter. So the actual observational possibility of 
the universality we found for astronomical objects should 
be further investigated. 
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